Black hole thermodynamics from simulations of lattice Yang-Mills theory 
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We report on lattice simulations of 16 supercharge SU (N) Yang-Mills quantum mechanics in the 
't Hooft limit. Maldacena duality conjectures that in this limit the theory is dual to IIA string 
theory, and in particular that the behavior of the thermal theory at low temperature is equivalent 
to that of certain black holes in IIA supergravity. Our simulations probe the low temperature regime 
for N < 5 and the intermediate and high temperature regimes for N < 12. We observe 't Hooft 
scaling and at low temperatures our results are consistent with the dual black hole prediction. The 
intermediate temperature range is dual to the Horowitz-Polchinski correspondence region, and our 
results are consistent with smooth behavior there. We include the Pfafnan phase arising from the 
fermions in our calculations where appropriate. 
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I. INTRODUCTION 

String theory has provided remarkable insight into the 
quantum physics underlying black holes. Much recent 
progress stems from conjectured dualities, which, in an 
appropriate limit, relate the finite temperature low en- 
ergy supergravity limit of the string theory to strongly 
coupled thermal field theory. The entropy of the black 
holes that arise in these supergravity theories can then 
be computed in principle by counting microstates in their 
dual field theories. The pioneering calculations of black 
hole entropy in [l], 0] are examples where the dual field 
theory is a 2-d conformal field theory which allows com- 
putation of the entropy despite the strong coupling. 

For a large number N of coincident Dp-branes in the 
'decoupling' limit [H, 0| the dual field theory is (1 + p)- 
dimensional strongly coupled maximally supersymmetric 
SU(N) Yang-Mills theory, taken in the 't Hooft limit. 
The case of D3-branes yields the AdS-CFT correspon- 
dence. Analytic calculation of the corresponding black 
hole entropy of these theories has proven elusive despite 
interesting attempts ||. 

Here we use lattice methods to study the thermal gauge 
theory and hence test these conjectured dualities. The 
simplest case for lattice work corresponds to DO-branes 
@, where the dual is thermal 16 supercharge Yang-Mills 
quantum mechanics (the 'BFSS model' 0]). This theory 
has recently been numerically studied using a non-lattice 
formulation 0, Q. Earlier analytic approaches used a 
variational method [10L I 111 . Related zero temperature 
numerical works are |l2tll3. [T^|. 

In this letter we simulate the super quantum mechan- 
ics in the t 'Hooft limit over a range of temperature and 
present preliminary results. We obtain intermediate tem- 
perature results for N < 12 and low temperature results 
for N < 5. We pay particular attention to the continuum 
limit and the behavior of the important Pfafnan phase 
arising from the fermions. More details of the method 
and results will be given in [l5j . 



II. DUALITY AND BLACK HOLES 

The type IIA string theory reduces to a supergrav- 
ity theory for low energies compared to the string scale 
(o/) -1 / 2 . In this limit the thermal theory contains black 
holes with N units of DO-charge. Their energy, E, is 
a function of their Hawking temperature, T. Defin- 
ing A = Ng s a'~ 3 ' 2 where g s is the string coupling, 
we may write a dimensionless energy and temperature 
e = EX^ 1 / 3 and t = TA~ 1,/3 . One finds provided we take 
N large and t « 1 the black hole is weakly curved on 
string scales and the quantum string corrections are sup- 
pressed. The energy of this black hole can be precisely 
computed by standard methods [4] giving, 
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Duality posits that the thermodynamics of this black hole 
should be reproduced by the dual Yang-Mills quantum 
mechanics at the same temperature with <? s o/~ 3 / 2 = g\ M 
so A is identified with the 't Hooft coupling. 

In the large N limit, at high temperatures t > 1, the 
bound state of DO-branes is of order the size of the string 
scale, and hence all a' corrections are important. One 
should best think of the configuration dominating the 
partition function as a hot gas of DO-branes bound by 
strings. Horowitz and Polchinski have argued that the 
low temperature black hole and high temperature gas 
are the asymptotic descriptions and intermediate tem- 
peratures smoothly interpolate between these fiH ]. 

III. LATTICE IMPLEMENTATION 

The 16 supercharge SU(N) Yang-Mills quantum me- 
chanics arises from dimensional reduction of J\f = 1 su- 
per Yang-Mills in 10-d. The 10-d gauge field reduces 
to the 1-d gauge field A and 9 scalars, X 1 , i = 1, ... ,9 
and the 10-d Majorana-Weyl fermion to 16 single com- 
ponent fermions, a = 1, . . . , 16. All fields transform 
in the adjoint of the gauge group. In order to simulate 
the theory we must integrate out the fermions giving rise 
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to a Pfaffian. The continuum Euclidean path integral, 



Z = J dAdXPf (O) < 



is then given by, 
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We take Euclidean time to have period R. 



The 7 T ,7* are the Euclidean Majorana-Weyl gamma 
matrices, and we choose a representation where, 7 1 " = 
Id 8 N 
Id 8 

We have a choice of fermion boundary conditions. 
Thermal boundary conditions correspond to taking the 
fcrmions antiperiodic on the Euclidean time circle and 
correspond to a temperature t = A -1 / 3 JR. We will also 
employ periodic fermions and the continuum partition 
function is then an index, with t the inverse volume. 
The Pfaffian is in general complex [17] . It is important 
in principle to include the phase of the Pfaffian in the 
Monte-Carlo simulation, and we discuss this later. 

We discrctizc this continuum model as, 
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where we have rescaled the fields X^ a and ^i, a by pow- 
ers of the lattice spacing a — R/L where L is the number 
of lattice points to render them dimensionless. We have 
introduced a Wilson gauge link field U a , and taken co- 
variant difference operators (D_W) a = W a —U%W a -\U a , 
(D + W) a = U a W a+1 W a - W a . Notice that the fermionic 
operator is free of doublers and is manifestly antisym- 
metric. This lattice action is finite in lattice perturba- 
tion theory and hence will flow without fine tuning to the 
correct supersymmetric continuum theory as the lattice 
spacing is reduced @, [l5| . 

We use the RHMC algorithm [l8[ to sample configura- 
tions using the absolute value of the Pfaffian. The phase 
may be re- incorporated in the expectation value of an ob- 



servable A by reweighting as < A >= 



E,„(-A^) 
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Here 



is the phase of the Pfaffian and the sum runs over 



all members of our phase quenched ensemble. 

We find in practice that the RHMC simulation of the 
thermal theory at low temperature, t < 1, exhibits an 
instability corresponding to the scalar fields moving out 
along the flat directions of the classical potential. Hence 
the algorithm never thermalizes and cannot be used to 
approximate the path integral. This has been observed 
before Q. We believe this divergence is a lattice arti- 
fact that is related to the discretization of the fermion 
operator. In previous work [6| we have simulated the 4 
supercharge quantum mechanics over a range of t, using 
a Weyl representation for the fermions where one ob- 
tains a real positive determinant. However, we have also 
tried using a Majorana representation where one obtains 



FIG. 1: Top: Plot showing e/t verses dimensionless temper- 
ature t for the periodic SU(2) theory for various numbers of 
lattice points. Bottom: Plot of the Polyakov loop against 
temperature for the same theory. 

a Pfaffian which we have discretized in analogy with the 
16 supercharge case discussed here. While no divergence 
of the scalars was observed in the Weyl simulations over a 
large range of t [f|, the Majorana implementation has the 
same instability we observe in the 16 supercharge case for 
t < 1. Since both representations are equivalent in the 
continuum limit this implies that the instability is not a 
property of the continuum theory as is claimed in Q but 
merely an artefact of finite lattice spacing. More details 
will be given in [la |. 

We find no such problem simulating the periodic the- 
ory at small t. At low temperature we expect the thermal 
and periodic theory to be similar, and the configurations 
that dominate the path integral will be similar. Hence in 
order to simulate the thermal theory at low temperature, 
t < 1, we have employed a reweighting of the periodic 
theory. We can expect to get good results for the ther- 
mal theory at low temperature by computing expectation 
values using the periodic theory, and reweighting as, 



< A > T = 



where Ylm^ is a sum over the phase quenched ensemble 
generated for the periodic theory, Op and Ot are the 
periodic and thermal fermion operators respectively, and 
< . . . >t is the expectation value for the thermal theory. 

IV. RESULTS 

We have simulated the thermal and periodic theories 
concentrating on the range 0.3 < t < 5. We have focused 
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FIG. 2: Top: Plot of dimensionless energy e/t verses dimen- 
sionless temperature t for the thermal SU(5) theory with 5 
and 10 lattice points. Bottom: Plot of the cosine of the Pfaf- 
fian argument for the thermal 517(5) theory with 5 points. 



on two observables - the mean energy, e, and absolute 
value of the trace of the Polyakov loop, P. In the Yang- 
Mills theory these are given by [f| [l~5| 



<e/t> = ^ (h(N 2 - 1)- < S bos > 
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p = ^<i^IKi> 



(5) 



a=0 



The inclusion of l/N 2 , 1/N in these definitions is to en- 
sure these quantities are finite in the t'Hooft limit for 
a deconfmed phase. In the periodic case, since Z is an 
index, it should not depend continuously on the inverse 
volume i, and hence in the continuum e = 0. 

To check for a restoration of supersymmetry we have 
computed e/t in the periodic theory for a variety of lat- 
tice sizes L = 5, 10 and 20. The upper plot of figure Q] 
shows e/t for SU(2). For large t the index e is already 
consistent with zero for L = 5, while at small t it appears 
to approach zero as L increases. Notice that while this 
quantity is a sensitive test of the restoration of supersym- 
metry in the continuum limit L — > oo other observables 
such as P shown in the lower plot are relatively insensi- 
tive to the number of lattice sites for L > 5. 

We have also examined the continuum limit of the ther- 
mal theory. In figure [5] we show L = 5 and 10 data for 
the thermal energy for SU (5) (in the phase quenched ap- 
proximation - which we discuss shortly). As noted above, 
we find a lattice instability for the thermal theory with 
t < 1 (with some dependence on N and L). However, 
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FIG. 3: Top: A plot of the dimensionless energy e/t verses 
dimensionless temperature t. Data shown is generated in two 
ways. For temperatures larger than t ~ 1 we simulate the 
thermal theory for N = 2, 3, 5, 8, 12 with 5 points. The low 
temperature results are computed for N = 2, 3, 5 for 5 points 
by simulating the periodic theory, and reweighting with the 
appropriate combination of the thermal and periodic Pfaffi- 
ans, as described in the main text. The low temperature black 
hole prediction is shown. Bottom: A plot of the Polyakov loop 
observable P for the same cases. 
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FIG. 4: For comparison with figure [3] e/t verses t is shown 
for the quenched theory for N = 5, 12, 30 and periodic theory 
with Pfafhan reweighting for iV = 3, 5, using 5 point lattices. 



for larger t this does not occur. As argued above we be- 
lieve this is an artifact of our lattice formulation and has 
nothing to do with continuum physics. The points plot- 
ted in the figure are taken only from simulations where 
the scalar distribution remained bounded for hundreds 
of physical RHMC times (the observed instability sets in 
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very quickly in RHMC time, so the change in behavior 
is easy to identify). The plot shows that these lattice 
spacing effects are small and hence for the remainder of 
our results we show only data from L — 5 point lattices. 

In the lower plot of figure [2] we show the mean cosine 
of the Pfaffian phase for the thermal SU (5) theory with 
L = 5 lattice sites. As expected this phase becomes more 
important at lower temperatures but the actual value is 
close to one over the range of temperatures where we can 
directly simulate the thermal theory. Indeed the effects 
of reweighting are negligible in this temperature regime. 
Hence for the data we present later for direct simulation 
of the thermal theory we use the phase quenched approx- 
imation. Since the Pfaffian is very costly to compute this 
allows us to work at larger N . 

We now turn to the main results of this letter. In 
figure [3] we plot the energy and Polyakov loop for various 
N and L = 5 point lattices. For high temperature we 
have used direct simulation of the thermal theory (phase 
quenched as discussed above), and we are able to obtain 
results up to TV = 12. At low temperatures we obtain 
results by a reweighting of the periodic simulations as 
discussed above. The results from both methods agree 
in the regime where they overlap t ~ 1. 

At very high temperatures the curves approach a con- 
stant corresponding to the result from classical equiparti- 
tion assuming ./V 2 deconfincd gluonic states (the fermions 
are lifted out of the dynamics by their thermal mass in 
this limit). In contrast for low temperatures the energy 
approaches zero signaling the presence of a supersymmet- 
ric vacuum at vanishing temperature. 

As seen before [f|, we see that t'Hooft scaling sets in 
for small N, with N = 3 already giving results close to 
an extrapolated large N result. The high temperature 
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asymptotics computed in [191 ] are also plotted and agree 
with our data. Our results also appear consistent with 
those found recently using non-lattice methods 

There are two important physical observations. Firstly 
the curves appear to interpolate from high to low tem- 
perature smoothly - there is no obviously discontinuous 
behavior. This is to be contrasted with the quenched 
version of this theory which has a large N confine- 
ment/deconfinement phase transition at t ~ 0.9 [2oL l2lj|. 
Since the intermediate temperature range t ~ 1 is dual 
to the regime where the thermal DO-branes have a ra- 
dius comparable to the string scale, we are probing the 
Horowitz-Polchinski correspondence regime, and seeing 
apparently smooth behavior there. 

Secondly, the low temperature behavior of the theory 
appears consistent with the prediction from supergravity, 
also shown in the plots. This is to be contrasted with the 
quenched energy curve shown for comparison in figure 
|4] which departs strongly from the black hole prediction 
at low temperature. In this figure we also show the pe- 
riodic theory which shows the degree of supersymmetry 
breaking for these lattices is small. 

It would be very interesting to extend these calcula- 
tions to 2 and 3 dimensional Yang-Mills systems which 
are thought to be dual to Dl and D2 brane systems using 
recent lattice formulations retaining exact supersymme- 
try [13. 
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